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Abstract 

In the final stages of collapse, quantum radiation due to particle creation from a 
naked singularity is expected to be significantly different from black hole radiation. 
In certain models of collapse it has been shown that, neglecting the back reaction 
of spacetime, the particle flux on future null infinity grows as the inverse square 
of the distance from the Cauchy horizon. This is to be contrasted with the flux of 
radiation from a black hole, which approaches a constant (inversely proportional 
to the square of its mass) in the neighborhood of its event horizon. The spectrum 
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of black hole radiation is identical to that of a black body at temperature T = 
(8irM)~ 1 . We derive the radiation spectrum for a naked singularity formed in 
the collapse of a marginally bound inhomogeneous dust cloud and show that the 
spectrum is not black body and admits no simple interpretation. 
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The singularity theorems of Penrose, Hawking and Geroch ensure that, un- 
der fairly general conditions, the collapse of a very massive star will end in the 
formation of a singularity The theorems do not by themselves indicate, however, 
whether the singularity will be covered by an event horizon or whether it will be vis- 
ible to the external observer. If the singularity is visible to the external observer, 
it is said to be naked. Naked singularities are further categorized according to 
whether or not they are visible to the asymptotic observer. Those that are visible 
asymptotically are globally naked, otherwise they are locally naked. Both varieties 
have long been considered undesirable for several reasons and this has given rise to 
the so-called "Cosmic Censorship Hypothesis" M (CCH) which essentially banishes 
them from the physical universe. While the general consensus seems to be in favor 
of the CCH, numerous studies of classical collapse in the literature show that, on 
the contrary, naked singularities do indeed form from reasonable initial data^ 2 ' in 
classical relativity and low energy string induced gravity, making it necessary to 
understand just how nature avoids them, if it does so at all. Indeed, it is possible 
that the very processes that nature employs to rid the universe of naked singu- 
larities may have astrophysically observable consequences, which, in turn, can be 
exploited to obtain experimental information on the behavior of quantum fields in 
extremely curved spacetimes and even of quantum gravity itself. 

If matter does indeed attempt to collapse into a shell-focusing naked singularity, 
we can expect its behavior to differ significantly from that of matter collapsing into 
a black hole. This is because, toward the end stages of the collapse, not only is 
a region of high curvature accessible to the external observer, making quantum 
effects such as particle production and, eventually, the back reaction of spacetime 
particularly important, but also because the causal structure of the spacetime is 
vastly different. As a consequence of the exposure of regions of very high curvature, 
a feature that appears both within the context of exactly solvable models of two 
dimensional dilaton gravity^ and of more traditional collapse models in Einstein 
gravityt 4 ' 5 ' 6 ^ is the divergence of the radiation flux from the singularity on the 
Cauchy horizon. This effect was first observed in two dimensional stress-tensor 
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calculations given in ref. [6]. On the contrary, black hole radiation approaches a 
steady state near the event horizon and shell-crossing naked singularities, which 
have been analyzed by Ford and Parker' 7 ], have been shown not to produce a large 
flux of particles. 

In earlier articles' 4 ' 5 ] , we had considered particle production due to the marginally| 
bound, self similar collapse of inhomogeneous dust in Einstein gravity. The model 
admits a globally naked singularity at the center when the mass parameter falls 
within a certain range and gives rise to a black hole otherwise. Our analysis estab- 
lished in two ways that the stress tensor of evaporation grows without limit as the 
Cauchy horizon is approached. Naturally, the evaporating cloud cannot radiate 
away more energy than it has, implying that, barring some unexpected effect due 
to the back reaction of spacetime, the Cauchy horizon can safely be expected not 
to form at all. Thus, the CCH (at least in its weak form) would have its origin 
in the quantum theory, not in the classical. Given that the radiation flux is large, 
particularly toward the end stages of the collapse, it is natural to expect that 
the effects of the flux may be observable. If so, what would be its signature? In 
search of the answer to this question, we will consider the spectrum of the radiation 
emitted from the collapsing cloud when the mass parameter is such that a naked 
singularity would have been classically formed. As it turns out, the spectrum is 
unique and has no simple interpretation in terms of known distributions. Contrary 
to the intensity of the radiation flux at infinity which, as we have said before, is 
due to the regions of high curvature encountered, the peculiar nature of the spec- 
trum we derive owes to the causal structure of the spacetime. It is particularly 
interesting as its very uniqueness may serve to distinguish experimentally between 
naked singularities and other celestial radiators. 

In the following, we have neglected the scattering of waves within the collaps- 
ing cloud and the back reaction of spacetime, confining ourselves to the geometric 
optics approximation. Clearly, both of these simplifications will imply that correc- 
tions to the spectrum obtained below are to be expected. These corrections will be 
frequency dependent and their importance will depend on how close the collapsing 
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cloud actually comes to forming the Cauchy horizon. The more important correc- 
tions to the spectrum will arise from the back reaction of the spacetime. In low 
energy string theory models of gravity, the time at which the back reaction becomes 
important is estimated by evaluating the size of the dilaton coupling constant at 
the center of the infalling matter. There is no such simple method in classical rel- 
ativity, but we might guess that it will become significant for wavelengths that are 
smaller than the smallest radius of curvature encountered though, unfortunately, 
without the benefit of a full quantum theory of gravity, it is impossible to estimate 
precisely the strength or nature of these corrections. 

Let us begin by briefly reviewing the classical collapse of a marginally bound 
inhomogeneous dust cloud. The matter is described by the stress energy tensor 

TS = e(t,r)5K (1) 

and the solution of Einstein's equations in comoving coordinates is the Tolman- 
Bondi^ metric 

ds 2 = dt 2 - R 2 {t,r)dr 2 - R{t,r) 2 dQ 2 (2) 

where 

/ I \ 2/3 

R(t,r) = r\l - (3) 

and R'(t, r) is the partial derivative of R(t, r) with respect to the coordinate r. The 
function F(r) is called the mass function. We will consider the special case of "self- 
similar" collapse which corresponds to a specific choice of F(r), namely F(r) = Xr 
where A is a positive constant and is related to the mass of the collapsing matter. 
The solution is obviously valid only in the interior of the cloud, that is up to a 
fixed value, r , of the coordinate r, giving the total mass as 2M = \r Q . Beyond 
r the metric is Schwarzschild and the two regions, the interior and the exterior, 
can be made to match smoothly at r = r . It can be shown that R(t, r) = is a 



5 



curvature singularity which gives the singularity curve in comoving coordinates as 
t = 2r/3 v / A. 

The solution above has been examined in detail in [4] and here we will only 
describe the spacetime. For self-similar collapse, the null coordinates inside the 
cloud are given as 



u 



+re 7 - x - R' > 
-re J ~ x - R! < 

(4) 

+re 1+ x + R' >0 
-re J + x + R' <0 

in terms of x — t/r and the integrals I± defined by 

They have been defined so as to reduce to the usual null coordinates of flat 
Minkowski spacetime in the limit as A — > 0. To analyze the causal structure it 
is convenient to use the variable y = yj R(t, r)/r in terms of which the integrals in 
(5) can be expressed as 

where f±(y) are the quartic polynomials 

f±(y) = 3y 4 =F ^y 3 — 3y + 3a/A (7) 

The center of the cloud can be shown to be just the line u = v , whereas the 
singularity curve away from the origin is spacelike and given in terms of the null 
coordinates of (4) by v — —cu where c is a positive constant. It is the origin, 
u = = v , that is of interest, the singularity being globally naked here if and only 
if there exists at least one positive real root of the polynomial f-(y) defined in 
(7). [4,8] Each positive real root of f-(y) corresponds to an out going null ray that 



6 



reaches X + , whereas each positive real root of the polynomial f+(y) corresponds 
to an infalling null ray originating on Z~ and intersecting the origin. Now it can 
be shown that f-(y) admits two positive real roots when the mass parameter is 
less that a critical value, A c ~ 0.18, and none when the mass parameter A > A c , 
while the polynomial f+(y) always admits a pair of real roots, one negative and 
one positive, when A is in this range. A c is thus a critical point. The singularity at 
the origin is globally naked for A < A c and for A > A c it is covered. In the latter 
case, the collapse leads to the formation of a black hole. 

The exterior region is described by the usual Schwarzschild metric 

-R{U,V)/k 

ds 2 = R{JJV) dUdV - R 2 (U,V)dQ 2 (8) 

where T, R are the Schwarzschild coordinates, Q is the solid angle and U, V are the 
Kruskal coordinates defined in the usual way by 

U = -2Ke-°/ 2K 

(9) 

V = 2ne v ' 2 \ 

in which definition k = 2M and U, V are the Eddington-Finkelstein outgoing and 
incoming null coordinates, U = T — R*,V = T + R* respectively. i£* is the tortoise 
coordinate. The two spacetimes in (8) and (2) are matched by matching the first 
and second fundamental forms at the boundary r = r J 7 l Comparing the angular 
parts of the metrics, it is natural to let R(t) = R(t, r ) at the boundary, whence 
one easily derives the relationship between the Schwarzschild null coordinates and 
the variable y on the boundary, 

U(y) = - ^=y 3 2V\r y - r y 2 - 2Ar ln \y/V\ - 1| 

Z X (10) 
V(y) = - ^=y 3 - 2V\r y + r y 2 + 2Ar ln \y/V\ + 1|. 

While it is difficult to obtain an analytical relation between the null coordinates in 
the interior and exterior regions all along the boundary, it is also not necessary. We 
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will be interested in relating the two only near the Cauchy horizon. The Penrose 
diagram for the classical formation of a naked singularity is shown in figure I below. 

r=2t/JJx (v= -cu) 
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Figure I: Formation of a naked singularity when A < A c 

We will consider a massless scalar field, <f), propagating in this background. 
In the infinite past, on X - , the field operator, can be expanded unambiguously 
in terms a complete, orthonormal family which contains only positive frequencies 
with respect to the canonical affine parameter on I" . This expansion defines the 
"in" vacuum. In the region close to where the Cauchy horizon would classically 
be formed, the massless field should also be determined by its data on the Cauchy 
horizon and on future null infinity, I + . Two problems are usually raised about 
this program, viz. (a) the absence of a well defined initial value problem to the 
future of the intersection of the Cauchy horizon with I + implies that there is no 
justification for assuming that the totality of I + exists, putting in question the 
completeness of the mode functions there and (b) because of the singularity at 
u — = v, the construction of a complete orthonormal basis set of infalling waves 
on the Cauchy horizon would be impossible. The former objection is less serious in 
that the spacetime to the future of the Cauchy horizon is the analytic continuation 
of the spacetime to its past and, in as much as no singularity actually intersects 
X + , there is no a priori impediment to constructing a complete basis set on all of 



X + . The latter objection appears more serious, but it is based on the assumption 
that the Cauchy horizon will actually form. A simple argument shows that it 
cannot. In refs. [4,5] the power radiated across future null infinity in this model 
was calculated and seen to grow, to leading order, as 



P(U) = 



48tt 



7 2 ~1 
{Uo-Uf 



(11) 



where U — U is the Cauchy horizon. There is every indication that the result 
is generic (i.e., independent of the actual collapse model considered as long as it 
leads classically to the formation of a naked singularity) and it was first suggested 
by two dimensional stress-tensor calculations in ref. [6]. Above, 7 is a parameter! 4 ! 
that depends in a complicated way on the mass parameter, A. It approaches unity 
when A — > and increases, approaching ~ 1.53 when A — > A c . The total radiated 
energy is the integrated flux and, because the collapsing cloud cannot radiate more 
energy than it possesses, an upper limit being provided by the total mass, M of 
the cloud, one may estimate the retarded time, U', at which the cloud will have 
given up all its energy to the Hawking radiation. It is 

& ~ °' ~ T^F < °" (12 > 

For a sufficiently large total mass, M, the collapsing cloud can approach arbitrarily 
close to forming the Cauchy horizon, but it will not actually do so. Note that we 
have neglected the back reaction in the estimate above. Naturally, this will become 
important toward the end stages of the evaporation, before the cloud has radiated 
away completely, and at this point the spacetime will begin to change significantly. 
This change will be such as to prevent the Cauchy horizon from forming. Therefore, 
for a complete, orthonormal set of solutions to the wave equation which are purely 
infalling in the semi-classical picture, it is more reasonable to consider a null surface 
parallel to and in the past of the classical Cauchy horizon. Such a choice would 
not be plagued by the difficulty in (b) above. Moreover, the precise choice of 
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surface and functions is irrelevant because the functions in this basis play the role 
of spectator and do not influence the final result. 

Let us imagine that a null ray U = const., when traced backwards is found to 
originate in the infalling null ray V = G(U). From the general theory, ^ one knows 
that the number distribution of Minkowski particles observed on T + is simply given 
by the Bogoliubov coefficient 

oo 

/3(u/,u>) = [ -JV^^gfjj) (13) 
J Attvlulu' 

— oo 

(where the integral is performed over all of X + ) as 

oo 

M {0\N(u)\0) M = J dw'\(3{u',u)\ 2 (14) 



where |0)m is the Minkowski vacuum. 

When a globally naked singularity is due to be formed in the classical theory, 
the putative Cauchy horizon lies in the retarded past of the event horizon at some 
finite value, U , of U and, as we mentioned before, no radiation can be expected 
beyond this point as the cloud will have already given up all of its energy. This 
is consistent with the semi-classical picture in which it is easy to see that no 
incoming rays on Z~ may turn into outgoing rays in the retarded future of the 
Cauchy horizon. As A < A c , f-(y) admits two positive (real) roots each of which 
corresponds to a null outgoing ray from the origin. Let a{ represent the real roots 
of f-(y). Because y = ai implies that t = 2r(l — a 3 )/3\/A, we take the larger of 
these two roots as the one that gives the earliest null ray from u = = v and call 
it a-. The Cauchy horizon is therefore given by (y = a_) u = in the interior 
and 

U = ^^a 3 - 2v / Ar a_ - r a 2 _ - 2Ar ln \oc-/V\ - 1| (15) 
2V A 

in the exterior. As before, consider a ray that leaves the origin near the Cauchy 
horizon and in its retarded past. Such a ray will intersect the boundary of the 
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y_ = y-a- = [-—) (17) 



cloud at some u = const, such that y is close to a_, say y = a_ + y_. For small 

y-, 

I-(y) ~ 7-lny- + (16) 
where 7_ = 3a£_/ f!_(ai-), the prime denoting a derivative w.r.t. y. This gives 

1/7- 

To) 

Up to linear order in y- one therefore finds that 

u ~ u + r_y_ = c/ + r_ f --J (is) 

where r_(a_) = — 2r a; 3 _/A/A(Q;- — \/A). For A < A c , T_ is negative and monoton- 
ically increasing as a function of A. Applying the same reasoning to the infalling 
ray, which in the interior is given by v — 0, one finds that any infalling ray close to 
v = and in its advanced past will intersect the boundary at y = a + + y + , where 
a + is the positive (real) root of the polynomial f+(y), such that 

> 1/74 



y + = y-a+ = [--) (19) 
where 7+ = 3a+/ f+(a+). Thus one has 



\ 1/7+ 

v ~ v - r+ I -^J (20) 



where r + (« + ) = — 2r Q a\/ \f\(a + + \/A) and V = V(ct+) (as given in (10). Com- 
bining (18) and (20) after "reflecting" at the center (u = v) gives the relationship 
between V and U as 

v = g(u) = a - b(u -uy u < u (21) 

where 7 = 7-/7+ and there are no incoming rays corresponding to outgoing rays 
in the retarded future of the Cauchy horizon. A — U and B(> 0) are irrelevant 
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constants that can be evaluated explicitly. When A = (there is no infalling 
matter), one finds that 7 = 1 = B so that G(U) = U as expected. Thus, the 
Bogoliubov coefficient, (3{uj',uj), when 7 > 1 is 



Uo 



e i"'B(Uo-ur (22) 



Changing variables to z — (U — U), one has 

00 

P(w',u)) = ^\f^i e ~ iU3 ' A ~ i0jt)o f dze iuJz e iuj ' Bz " (23) 



which gives 

k=0 

The radiation from the naked singularity is clearly not black body radiation, as it 
is for a black hole. Indeed we are unable to find any analogy between (31) and the 
standard distributions that arise in particle physics. 

The expression above is useful to analyze the high frequency limit (Su'(w)" 7 — > 
0) limit of the spectrum, for in this limit it is sufficient to consider only the first term 
in the series. Integration over u' then yields the familiar logarithmic divergence in 
the high frequency region. The alternative expression 



{u',u)\ 2 



^ {iu>{jjB)-ihe-™ l^) k k + 1 2 
4^7^(^5)2/7 k\ 1 7 



serves to analyze its low frequency (Buj'iu))^ 1 — > 00) behavior. Integration over 
uj 1 in this limit shows a power law divergence in the infrared. This divergence is 
associated with the fact that there are an infinite number of quanta in each mode 
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on X + . As the collapsing ball produces a steady flux of radiation to Z + , the net 
flux for all time is infinite. The difference between the divergence in the low and 
high frequency regimes may be associated with the red-shifting of modes in the 
proximity of the putative Cauchy horizon. Nevertheless, \/3(u', to) | 2 is seen to be 
well behaved as a function of u. It approaches uj in the infrared and 1/cu in the 
ultraviolet regions. 

Both the low frequency behavior of the black hole and the naked singularity 
as well as their high frequency behavior are seen to differ significantly. At low 
frequencies the intensity of the radiation from a naked singularity drops off as u, 
whereas for a black hole it behaves as 1/cu. In the same way, at high frequencies, 
the radiation spectrum of a naked singularity drops off in intensity as 1/cu whereas 
the black hole spectrum drops off much faster, as e~^ HUJ where (3h is the inverse 
Hawking temperature. Naked singularities prefer to radiate at higher frequencies. 

In this paper we have described a potentially important astrophysical problem, 
namely that of a gravitational collapse that classically leads to the formation of a 
naked singularity. We have argued that the presence of strong gravitational fields 
accessible to the asymptotic observer imply that the naked singularity is likely not 
to form because the matter energy of the collapsing cloud will disperse as Hawking 
radiation before the arrival of the Cauchy horizon. This radiation will be intense 
and possibly observable asymptotically. We have also shown that the unique causal 
structure associated with a global naked singularity leads to a unique spectrum. 

In our treatment, we have made no attempt to include the back reaction of 
the spacetime though it may be argued that it is precisely the back reaction that 
will dominate the evolution of the system in the final stages of collapse. The back 
reaction is a serious issue and should be treated satisfactorily at some point. As 
far as we know, only string theory or its generalizations provide us with consistent 
quantum theories of gravity, so it would appear that it can be satisfactorily studied 
only in the context of string gravity. Yet, in the few known exactly solvable models 
of collapse in string theory, I 3 ' 13 ! the findings have not so far contradicted the general 



13 



conclusions which have been arrived at by the semi-classical approach. It is possible 
therefore that the above description is reasonably accurate. 

The subject of gravitational collapse in Einstein gravity is rich and cannot be 
ignored. The purpose of this work is to initiate a discussion on the subject and to 
attempt to find observational signatures that distinguish between matter collapsing 
to form a black hole and matter collapsing to form a naked singularity, even though 
that end state may never be achieved. (Black hole studies, for instance, are very 
useful even though the end state of the black hole is not known.) Our hope is that 
whatever information one can extract from models of collapse in Einstein gravity 
may later be incorporated into a better understanding of the role of the spacetime 
from string (or, for example, M— ) theory. 

To conclude, we venture to summarize in terms of a Penrose diagram (figure II) 
what we believe the spacetime in the fully quantized theory should look like, based 
upon the understanding that has been obtained from the semi-classical findings 
related above. 
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Figure II: Conjectured structure of the spacetime in the full quan- 
tum theory of a cloud attempting to form a classical naked singularity 

In Figure II, region I is described by the Tolman-Bondi metric given in eq.(2). 
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Region II is the Schwarzschild exterior and, in region III, Hawking radiation be- 
comes significant so that energy streams out of the cloud toward infinity. Strictly 
speaking, region III also contains collapsing matter, but the radiation dominates. 
As such it is probably acceptable to model this patch by the Vaidya metric. It 
must, of course, be matched to regions I and II in the usual way. It is along X + in 
this portion of spacetime that the spectrum obtained above should be a good ap- 
proximation until the end stages when the back reaction kicks in. Notice that the 
putative Cauchy horizon is to the future of this region. It falls in a patch (region 
IV) that is Minkowskian because, presumably, all the energy has by now escaped 
to infinity. This scenario will be discussed and the limitations of the approximation 
above will be systematically treated more fully elsewhere. 
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